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We present a measurement of the double longitudinal spin asymmetry in inclusive ir° production 
in polarized proton-proton collisions at y/s — 200 GeV. The data were taken at the Relativistic 
Heavy Ion Collider with average beam polarizations of 27%. The measurements are the first of a 
program to study the longitudinal spin structure of the proton, using strongly interacting probes, at 
collider energies. The asymmetry is presented for transverse momenta 1-5 GeV/c at mid- rapidity, 
where next-to-leading order perturbative quantum chromodynamic (NLO pQCD) calculations well 
describe the unpolarized cross section. The observed asymmetry is small and is compared with a 
NLO pQCD calculation with a range of polarized gluon distributions. 
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From polarized lepton-nucleon deep inelastic scatter- 
ing (DIS) experiments over the past 20 years it is known 
that only ~ 25% of the proton spin is carried by quarks 
and anti-quarks yj. The rest of the proton spin must 
hence be carried by the gluons and orbital angular mo- 
mentum. DIS experiments have constrained the possible 
gluon polarization in the proton through the measure- 
ment of scaling violation in inclusive polarized scatter- 
ing 2] , and through semi- inclusive measurements of two 
hadrons to utilize the photon-gluon fusion process Q. 
A fixed target experiment at Fermilab first presented a 
measurement with strongly interacting probes Q. The 
reach of these measurements was limited, due to the low 
energy available for fixed target experiments. Presently, 
the gluon contribution to the proton spin is largely un- 
known. 

The polarized proton collisions at the Relativistic 
Heavy Ion Collider (RHIC) provide a new laboratory to 
study the proton spin structure with strongly interacting 
probes. The PHENIX experiment has reported the un- 
polarized cross section for ir° production at mid-rapidity 
for pt—I-U GeV/c, which is described well by next-to- 
leading-order perturbative QCD (NLO pQCD) calcula- 
tions over 8 orders of magnitude |5j . In this Letter we re- 
port the first results on the double spin asymmetry All, 
for inclusive 7r° production at mid-rapidity in longitudi- 
nally polarized proton-proton collisions corresponding to 
0.22 pb _1 integrated luminosity with the PHENIX de- 
tector. 

In perturbative QCD All is directly sensitive to 
the polarized gluon distribution function in the proton 
through gluon-gluon and gluon-quark subprocesses 0. 

The double spin asymmetry in 7r° production is given 

by 



(1) 



where a ++ (<j-| ) is the cross section of the reaction when 

two colliding particles have the same (opposite) helicity. 
Here we neglect the parity violating difference in cross 

section between (++) «-* ( ) and (H — ) <-> ( — h) beam 

helicity configurations. Since the cross section can be 
obtained by dividing the experimental yield (N) by the 
integrated luminosity (L), All is expressed as 
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where Py(B) are the polarizations of the RHIC "yellow" 
("blue") beams, and R is the ratio of luminosities of pro- 
tons colliding with like to unlike hclicities. 

For the 2002-2003 RHIC run, 55 bunches of polarized 
protons, typically 5 x 10 10 protons per bunch, were loaded 
into each of the yellow and blue counter-rotating acceler- 
ator/storage rings of RHIC and accelerated to 100 GeV. 
The bunch lengths and separations were ~1 ns and 213 



ns, respectively. The beam polarization sign for each 
bunch was prepared independently at the source, with 
the successive bunches in one ring alternating in po- 
larization sign, and with successive pairs of bunches in 
the other ring alternating in sign. The locations of the 
bunches were identified relative to a RHIC timing clock. 
In this way, the experiments collected data from collisions 
with all four combinations of blue-yellow ring beam po- 
larization signs simultaneously. 

The stable direction of the proton spin in RHIC is ver- 
tical, but the region around the PHENIX experiment in- 
cludes sets of magnets (spin rotators) to rotate the spin to 
the longitudinal direction at the collision point, and then 
back to vertical after the interaction point, in order to 
provide collisions with longitudinal polarization, and to 
maintain the required vertical polarization around RHIC. 
The RHIC polarimeters measure the transverse beam po- 
larization away from the interaction points, independent 
of the operation of the spin rotators. 

The transverse beam polarization was measured in 
RHIC independently in each beam using proton-carbon 
elastic scattering in the coulomb nuclear interference 
(CNI) region (tJ. The analyzing power A p ^ was mea- 
sured for 22 GeV beam energy, A^ c (22), to ±30% 0. 
The energy dependence of the analyzing power over the 
RHIC energies is expected to be small, <10% [&]. For 
the results reported here, we have used the same an- 
alyzing power at 100 GeV as at 22 GeV, and 10% is 
added in quadrature to 8A P ^ (22) to give a ±32% uncer- 
tainty for ^4^(100). With these assumptions, the aver- 
age polarization in the analyzed data set in this paper was 
^/(P B P Y ) = [27.0 ± 0.3(stat) ± 3(syst) ± 8(A P ° syst)]%. 

Local polarimeters, sensitive to the transverse polar- 
ization at collision, were used to set up the spin rotators, 
and to monitor the beam polarization direction at the 
PHENIX experiment. The local polarimeters utilized a 
transverse single spin asymmetry in neutron production 
in p-p collisions at y/s = 200 GeV 10]. For vertically 
polarized beam a left-right asymmetry is observed for 
neutrons produced at very forward angles, with no asym- 
metry for production at very backward angles. A fully 
longitudinally polarized beam produces no asymmetry. 

Neutrons with E n > 20 GeV and production angle 
0.3 < 8 n < 2.5 mrad were observed by two hadronic 
calorimeters located ±18 meters from the interaction 
point (ZDC or Zero Degree Calorimeter |ll|) . Scin- 
tillator hodoscopes at 1.7 interaction length provided 
the neutron position at the ZDC, and thus the neutron 
production angle and azimuthal angle <p — arctan(a;/y) 
with y vertically upward. The x axis forms a right 
handed coordinate system with the z axis defined by the 
beam direction for forward production. The single-spin 
asymmetry e was calculated versus azimuth, from the 
four rates iV T) 0, i\^^ +7r ,iV|^,iV|^ +7r , using the geomet- 
ric mean |l2| . This method largely cancels differences in 
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FIG. 1: Single spin raw asymmetry normalized by the beam 
polarization, e/P, as a function of azimuthal angle <f>, for for- 
ward neutron production. The solid points and the curve 
correspond to the spin rotators off (transverse polarization) 
condition and the open points and dashed curve correspond 
to the spin rotators on (longitudinal polarization). Curves are 
sine function fits to the data, representing possible transverse 
polarization. The data are for individual runs, where the 
blue (yellow) polarization was 24% and 33% (8% and 28%), 
for spin rotators off and on, correspondingly. 



luminosity between ] and [ polarization collisions and 
between detector acceptance differences at <ft and + tt. 
Fig. 1 shows the observed asymmetry, for the spin rota- 
tors off and on, for the blue and yellow beams. With the 
spin rotators off, a left-right asymmetry is observed from 
the vertically polarized beam. With the spin rotators on, 
the measured transverse polarization, averaged over the 
run, was (P Bx ) = 3.3% ± 1.9%, (P By ) = 0.8% ± 2.0%, 
(P Yx ) = -2.0% ± 1.3% and (P Yy ) = 5.4% ± 1.7%, out 
of (P) = 27%. The double spin transverse polariza- 
tion was (P Bx Pyx) = (0.4 ±1.1) • 10~ 3 and (P By PYy) = 
(-0.2 ± 0.8) • 10~ 3 , compared to (P B P Y ) = 0.07. There- 
fore, with the spin rotators on, the transverse asymmetry 
is greatly reduced, indicating a high degree of longitudi- 
nal polarization. 

A separate run with the spin rotators set to give radial 
polarization confirmed the direction of the polarization 
for each beam. 

Collisions in PHENIX are defined by the coincidence 
of signals in two beam-beam counters (BBC) [J 3j located 
at ±1.44 m from the nominal interaction point and sub- 
tending a pseudo-rapidity range ±(3.0 — 3.9) with full 
azimuthal coverage. The BBCs select about half of the 
inelastic proton-proton collisions || . The vertex was re- 
constructed from the time difference of the hits in the 
two BBCs. The collision vertex was required to be within 
30 cm of the nominal interaction point. Events satisfying 
this condition constitute the minimum bias (MB) trigger, 
which was used for relative luminosity measurements. 

A coincidence of the two ZDCs was used to estimate 
the possible bias in the relative luminosity measurement 
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FIG. 2: Two photon invariant mass, M 77 , distribution for 
the pr = 2-3 GeV/c range. The n° signal was selected for 
±25 MeV from the it mass, indicated in black. The grey 
sidebands were used to estimate the All of the background 
under the tt° peak. 



from the BBCs. This was done by comparing the ac- 
cumulated number of triggers in the ZDCs and BBCs 
for each bunch and each fill. The accuracy of relative 
luminosity measurements 5R (Eq. 2) was estimated to 
be 2.5 x 10 -4 , which for the average beam polarization 
of 27%, translated to 5A LL = 1-8 X 10~ 3 , and, on the 
same uncertainty level, confirmed no A LL asymmetry of 
BBC triggers relative to ZDC. The ratio R averaged over 
the data sample used in the analysis was within 0.5% of 
unity. 



Neutral pions were reconstructed from the tt 
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decays using finely granulated (A(t> x A77 ~ 0.01x0.01) 
electromagnetic calorimeters (EMCal) [l4|, which con- 
sisted of two subsystems: a lead scintillator (PbSc) and 
a lead glass (PbGl) calorimeter covering three quarters 
and one quarter of the EMCal acceptance, respectively. 
Located at a radial distance of ^5 m from the beam line, 
the EMCal covered the pseudorapidity range of < 0.35 
and two azimuthal angle intervals of A0«9O° separated 
by <j> w 70° (nearly back-to-back). 

High pr 7r°'s were collected using coincidences between 
a MB trigger and an EMCal-based high pr photon trig- 
ger 0| . The trigger efficiency for 7r°'s varied from 8% in 
the 1-2 GeV/c p T bin to 90% in the 4-5 GeV/c p T bin. 

The 7T° reconstruction and photon identification cuts 
were optimized to minimize the background contribu- 
tion under the ir° peak in the invariant mass distribu- 
tion while keeping the ir° efficiency high. For photon 
identification we used the shower profile and the time- 
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of-flight measured by the EMCal, and charge veto cuts. 
The charge veto was set for those EMCal clusters asso- 
ciated with hit(s) in the pad chamber |l5j, which was 
located ~ 20 cm in front of the EMCal surface. In order 
to avoid the effects of electronic noise and suppress the 
very low energy background, only clusters with energy 
greater than 0.1 GeV in PbSc and 0.2 GeV in PbGl were 
used in the analysis. 

The ir° yield was extracted by integrating the two pho- 
ton invariant mass spectrum over a ±25 MeV/c 2 region 
around the ir° mass (signal region) as shown in Fig. 2 
by the dark band. The EMCal resolution was such that 
the widths of the 7r° mass peaks varied from 12 MeV/c 2 
in 1-2 GeV/c p T bin to 9.5 MeV/c 2 in 4-5 GeV/c p T 
bin, in both PbSc and PbGl. In the pr range of 1 to 
5 GeV/c, 4 million 7r° candidates were collected. The 
background contribution (combinatorial+hadronic) un- 
der the 7T° peak, r, varied from 27% in the 1-2 GeV/c 
bin to 8% in the 4-5 GeV/c bin, as indicated in Tabled 
The 7T° reconstruction efficiency due to photon identifi- 
cation cuts varied from 84% in the lowest pt bin to 93% 
in the highest px bin. 

The asymmetry of the background in the signal region, 
All > was evaluated using the asymmetry calculated from 
the grey bands on both sides of the signal (Fig. 2). The 
measured it asymmetry, A™™ j was corrected for the con- 
tribution of background using 
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The spin asymmetry for each beam fill |l6j A^ l L l was 
calculated using Eq.H For the A f ™ error evaluation, we 

considered only the N ++ and iV_| statistical errors. The 

resulting All was obtained after fitting a constant to all 
Aj^'s. The fit xju an d a "bunch shuffling" technique 
were used to check the uncertainties assigned to All. In 
each "bunch shuffling" we randomly assigned the helicity 
sign to every bunch crossing, keeping the balance between 
the number of bunches with correctly and inversely as- 
signed helicitics, so that the average polarization for each 
shuffled sample was nearly zero, and recalculated All- 
The widths of the distributions of All values obtained 
in all bunch shuffles were consistent with errors assigned 
to All indicating that all non-correlated bunch-to- bunch 
and fill-to-fill systematic errors were much smaller than 
the 7T° yield statistical errors. 

A number of systematic checks, including variation of 
photon identification criteria and mass window range for 
7r°s and background, were performed to look for possible 
systematic effects on the measured All values. None 
were found. 

The double spin asymmetries between (++) and ( ) 

and between (-1 — ) and ( — h) helicity configurations, as 
well as the single spin asymmetries for each polarized 
beam (Al — — IJ +~ cr - ) were evaluated. These measure 



TABLE I: Double spin asymmetry for the raw signal (it ± 
BG), for the background (BG) and for 7r° background cor- 
rected; single spin asymmetry for 7r° background corrected; 
for the four pr bins, the mean pr of 7r°'s being 1.59, 2.39, 
3.37 and 4.38 GeV/c, in four pt bins correspondingly . The 
numbers in brackets in the second column indicate the back- 
ground contribution to the 7T° Si gnal. 
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-1.5 ±0.9 (27%) 


1.6 ±1.4 


-2.7 ± 1.3 


-0.2 ±0.3 


2-3 


-1.5 ±1.1 (15%) 


-3.0 ±2.4 


-1.3 ±1.3 


-0.1 ±0.3 
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-1.8 ±2.5 (9%) 


-2.4 ±6.8 
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24 ± 17 
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FIG. 3: The measured double spin asymmetry A^ L versus 
mean pr of 7r°'s in each bin. A scale uncertainty of ±65% 
is not included. Two theoretical calculations based on NLO 
pQCD are also shown for comparison with the data (see text 
for details). 



parity violating asymmetries, if any. All of these asym- 
metries were consistent with zero. 

The results of the asymmetry are presented in Table [|] 
and Fig. [3] Systematic uncertainties for the asymmetry 
measurement are negligible. A total scale uncertainty of 
±65%, from the correlated polarization analyzing power 
uncertainty, 5A P ^ , for the two beams and the uncorre- 
cted measurement uncertainties, is not shown. 

Two theoretical curves based on NLO pQCD are shown 
in Fig. |3 representing different assumptions for the gluon 
polarization, with GRSV-std using the best global fit to 
inclusive DIS data, and GRSV-max using a gluon polar- 
ized distribution equal to the unpolarized distribution at 
the input scale of Q 2 = 0.6 GeV 2 [HOH. The ran g e 01 



theory curves reflects the uncertainty on gluon polariza- 
tion from the inclusive DIS measurements [lil H^ . |20| . 
The uncertainties from this measurement are similar in 
size to this range j^. Both curves show positive All- It 
was recently argued that a negative All would be diffi- 
cult to accommodate [i^ . The present results are consis- 
tent with the GRSV-std curve, with CL=16-20%, for the 
range in polarization scale uncertainty. For GRSV-max, 
CL=0.02-5%. These confidence levels do not include a 
theoretical uncertainty. In addition, the lowest px data 
point, at (pt)=1-59 GeV/c, may have a significant soft- 
physics contribution. The agreement of the NLO pQCD 
calculation for the cross section can only be checked 
to within the uncertainties of the calculation, estimated 
to be a factor two at this pr 

In summary, the reported results of the double helic- 
ity asymmetries for ir° production begin to explore the 
gluon polarization in the proton, using strongly interact- 
ing probes. The observed asymmetry is small, and the 
level of uncertainty of these data is comparable to that 
of the theoretical calculations adopting different sets of 
gluon polarizations deduced from global fits to the exist- 
ing deep inelastic scattering data. 
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